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Forced Oscillations in Fluid Tori and Quasi-Periodic Oscillations 
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Abstract. The kilo-Hertz Quasi-Periodic Oscillations in X-ray binaries could originate within the accretion flow, and be 
a signature of non-linear fluid oscillations and mode coupling in strong gravity. The possibility to decipher these systems 
will impact our knowledge of fundamental parameters such as the neutron star mass, radius, and spin. Thus they offer the 
possibility to constrain the nuclear equation of state and the rotation parameter of stellar-mass black holes. We review the 
general properties of these oscillations from a hydrodynamical point of view, when the accretion flow is subject to external 
perturbations and summarize recent results. 
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1. Introduction 

The modulations of the X-ray flux in Low-Mass X-ray Bina- 
ries (LMXBs) and their observed range in frequencies, from 
« 300 Hz to « 1200 Hz (see van der Khs 2000; McClin- 
tock & Remillard 2006; van der Klis, these proceedings, for 
reviews) strongly suggest that dynamical time scales in the 
inner regions of the accretion disks occurring in these sys- 
tems are being probed (i.e., at orbits only a few gravitational 
radii in circumference). 

In Black Hole X-ray Binaries (BHXRBs), the high fre- 
quency Quasi-Periodic Oscillations (QPOs) in the X-ray flux 
seem to be fairly stable, exhibiting nearly constant frequen- 
cies. Additionally, in the four galactic microquasars know to 
show more than one peak, their frequencies are in a 3:2 ratio, 
strongly suggesting that a resonance is responsible for their 
production (Kluzniak & Abramowicz 2000). 

In Neutron Star (NS) systems, the twin peaks observed 
in the frequency spectrum drift over a considerable inter- 
val (a few hundred Hz, as mentioned above), yet are lin- 
early correlated throughout this range. In many sources the 
peak separation, while not constant, is consistent with be- 
ing half the burst oscillation frequency, and in other sources 
with the burst frequency directly. In the one observed instance 
where burst oscillations, the spin frequency of the pulsar, 
and twin QPO peaks in the kHz range have been observed 
in the same source (SAX J1808.4-3658), the burst oscifla- 
tion is the spin frequency, and is twice the peak separation. 
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i.e., Pburst = Vs = 2Ai/(QPO) (Wijnands et al. 2003). In 
XTEJ 1807-294, the peak separation is consistent with the 
spin frequency (Linares et al. 2005), and there generally ap- 
pears to be a trend whereby the "slow" neutron stars have 
twin peaks with Ai^ w Vs and the "fast" ones show Av k, i/g 
(with the transition at about 300 Hz). 

We have previously suggested (Kluzniak et al. 2004; Lee, 
Abramowicz & Kluzniak 2004) that the peculiar behavior of 
SAX J 1808 can be understood in terms of a nonlinear re- 
sponse of the accretion disk, when subject to an external per- 
turbation at a fixed frequency. While it is currently not possi- 
ble to study the detailed structure and dynamical behavior of 
the accretion disk and its modes of oscillations in full detail, 
general physical principles can still be applied, and yield sug- 
gestive results. We have proposed that within the accretion 
disk, as in other nonlinear systems, a subharmonic response 
will appear in the presence of a perturbation, as will higher 
harmonics of the perturbation itself. Specifically, a second 
QPO frequency is to appear when two possible modes of os- 
cillation are separated by half the perturbing frequency, and 
thus couple. 

This presentation is devoted to a numerical study of 
the non-linear oscillations in such systems, by considering 
bound fluid tori representing a portion of the accretion flow 
around the compact object (in particular thin tori). We believe 
strong gravity plays an important role in their behavior, and 
consider its effects in a simplified approach. This is a purely 
hydrodynamical study, and so the potentially important ef- 
fects of the magnetic field have not been considered for now. 
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2.2 Fluid equilibrium configurations 



2. Hydrostatic equilibrium for fluid tori 



Astrophysical fluid bodies often appear in two topologicaUy 
different varieties: spheres and tori. In the former, support 
against gravitational collapse is mostly provided by pressure 
from within (e.g., by photons, degenerate or ideal gas par- 
ticles, neutrinos), with a possible centrifugal contribution, 
while in the latter, it comes mostly from rotation, with pres- 
sure providing corrections from simple test particle motion. 
Each of these terms: hydrodynamical and mechanical in na- 
ture, play an important role in the behavior of the system un- 
der external perturbations, and on the propagation of waves 
within such systems. As we will consider here the behavior 
and possible observational consequences of accretion disks 
in the particular context of LMXBs, we will focus mainly 
on toroidal configurations. Nevertheless, knowledge gathered 
from quasi-spherical systems is relevant, and indeed quite 
useful for interpreting these systems. 

2.1. Gravitational potentials 

At this point we may consider various forms for the grav- 
itational potential of the central mass, M. In the Newto- 
nian regime, obviously $ = <i>N = —GM/r. Since we 
are interested in applications to systems containing com- 
pact objects such as Neutron Stars and Black Holes, it is 
useful to consider pseudo-Newtonian potentials of the kind 
proposed originally by Paczynski and Wiita (1980), namely 
$pw = GM/{r - Tg), where Vg = 2GM/c^ is the grav- 
itational, or Schwarszchild, radius of the central mass. The 
important feature of this class of potentials is that they re- 
produce the existence of a marginally stable orbit, and that 
of capture orbits for finite values of the orbital angular mo- 
mentum, I, both unique features of General Relativity. In ad- 
dition to the original form, we have also used extensively a 
new pseudo-potential, $kl = GM[1 — exp(r„is/r)]/r„is, 
where r„is — ivg (Kluzniak & Lee 2002). The main advan- 
tage of this expression is related to the epicyclic frequencies 
for test particle motion. 

In Newtonian gravity, a test particle in circular motion 
with frequency around a point mass, M, will perform 
small epicyclic motions when perturbed. These can be ra- 
dial or vertical, but will always occur at frequencies Vr — 
Vz = — Q,/2n. This is just a re-statement of the well- 
known fact that the Newtonian orbits of a point mass are 
closed curves. In addition, no matter what the value of an- 
gular momentum a particle has, I, one can always place it a a 
certain radius, rcirc — P' jGM such that it will be in circular 
orbit. 

In strong gravity this is no longer the case, as the cen- 
tral well in the potential is so powerful that two qualitatively 
new effects appear. First, capture orbits exist even at finite I. 
Second, in the most simple case of the Schwarzschild met- 
ric (static, spherically symmetric geometry), the three-fold 
degeneracy between orbital and epicyclic frequencies is bro- 
ken, such that Vr <Vz — v^, and i/^ < for r < r,ns. Radial 
oscillations are thus unstable inside r,ns, and no circular orbit 
is possible in that region. 




Fig. 1. Effective potential wells for the Newtonian (N), Pa- 
cyznski- Wiita (PW) and Kluzniak-Lee (KL) formulae given 
in the text, in units of The value of angular momen- 
tum is given in units of VgC. Capture orbits are not possible 
in Newtonian physics, but they appear qualitatively in the 
pseudo-Newtonian potentials, as in General Relativity. The 
local minimum in each curve gives the radius of the circular 
orbit possible at the corresponding value of Iq. 



The particular form of $kl is such that the ratio of ra- 
dial to vertical epicyclic frequencies as a function of radius, 
Vrjvz is exactly equal to the value computed in full General 
Relativity in the Schwarszchild metric. Specifically, we have 

'~^*^exp(rins/r), (1) 
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Figure n shows the effective potential well as a function 
of radius for the three forms mentioned here and the same 
value of the orbital angular momentum. 

2.2. Fluid equilibrium configurations 

The fluid constituting the torus is located in the vicinity of the 
mass M, which produces the potential For simplicity we 
shall neglect the contribution of the mass density of the fluid, 
p, to the potential, and assume azimuthal symmetry, so that 
the problem can be studied in cyUndrical coordinates (r, z). 
The gravitational pull from the mass M is countered by the 
fluid pressure gradients and the centrifugal force, such that 

- -V$eff, (3) 
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3 NUMERICAL METHOD AND INITIAL CONDITIONS 
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Fig. 2. Meridional cross section of torus boundaries for a 
constant distribution of angular momentum, £ and differ- 
ent values of $o in the Paczynski-Wiita potential. In the 
limit of a slender torus (top) the meridional cross sections 
of the equipotential surfaces become ellipses with axis ratio 
flr/oz = Vr/vz- In the limit of a thick torus (bottom), the 
outer surface becomes spherical, while a narrow funnel de- 
velops along the rotation axis. 



is the effective potential and I is the distribution of specific 
angular momentum of the fluid, which depends only on the 
radial coordinate (according to Von Zeipel's theorem, for a 
polytropic equation of state the surfaces of constant angular 
velocity for a rotating fluid are cylinders). 

Now, equation (|3} can be trivially integrated over r if the 
fluid obeys a polytropic relation of the kind P — Kp^ , with 
K a constant and 7 the adiabatic index, to give 
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The constant of integration $0 is related to the size of the 
torus. The boundary of the torus is the surface over which the 
pressure is zero, and it coincides with a particular equipoten- 
tial surface (see Figure|2|l. 

The fluid within a torus is in hydrostatic equilibrium, with 
the pressure gradient balancing gravity vertically, and cen- 
trifugal forces providing the bulk of the radial support against 
the gravitational field. The circle of maximum density (or 
pressure) defines the center of the torus at rg. At this point. 



the pressure gradient vanishes, and thus the rotation is exactly 
Keplerian, as for a free particle. If the potential well is filled 
only slightly with gas, the torus will be slender, and confined 
to regions close to the equatorial plane. If the potential well 
is increasingly filled, the torus will grow and deform until it 
overflows the inner Roche lobe and accretion onto the cen- 
tral mass begins. This Roche lobe is analogous to that which 
occurs in binary systems, except in this case it is an effect 
of General Relativity (Abramowicz, Calvani & Nobili 1983), 
specifically of the existence of capture orbits for a finite value 
of angular momentum (see Figure^. 

It is easiest to think of the rotation law through the distri- 
bution of angular momentum, and convenient to write it as a 
power law, with £(r) = ^o('"/''o)"- For Keplerian rotation in 
a Newtonian potential we have, obviously, a = 1/2, while a 
constant distribution of angular momentum, such as the one 
used to plot the curves in Figures ^ and |2l has a = 0. For 
sub-Keplerian rotation (i.e., a < 1/2 in this case), the torus 
is bound by a surface of zero pressure, becoming infinite in 
the exact Keplerian limit (where one would have pressureless 
dust in orbit). The overall shape and size of the torus is thus 
determined by the degree to which the potential welU is filled, 
and by the rotation law. 

3. Numerical Method and Initial Conditions 
3.1. Numerical Method 

We have used the SPH (Smooth Particle Hydrodynamics) 
numerical method to perform our simulations (Monaghan 
1992). This is a Lagrangian formulation, which interpolates 
a given function A{r), with a set of given quantities A{r') 
through a kernel function W{r, h), using the convolution in- 
tegral: 



A{r) = / A{r')W{r - r')dr. 



(6) 



We have implemented the prescriptions of Monaghan & Lat- 
tanzio (1985) for azimuthal symmetry, using the kernel: 

3(7r)V2 + 3a)V4 ifO<|,<l; 
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Here h represents the smoothing length, and is comparable to 
the typical separation between fluid elements (it essentially 
gives the spatial resolution of the calculation), and r is the 
radial coordinate. In two dimensions, v — 2 and a — 1Q/7tt. 
The gas is modeled as an inviscid fluid, and so the Navier- 
Stokes equations take the form: 
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where R = Vr^-T^ is the distance to the central mass 
M. The sub-index art indicates the artificial viscosity terms. 



4 FORCED OSCILLATIONS 



which is used to account for the presence of shocks and to 
avoid interpenetration of the fluid elements. 
The energy equation is: 

du /P' 



dt 
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where u is the internal energy per unit mass. No external (i.e., 
radiative) losses are assumed, and thus the only change in u 
arises from work done. When discretized over a finite number 
of fluid elements, often called particles in SPH, the convolu- 
tion integral becomes a sum over all elements. 

We have used the prescription given by Balsara (1995) for 
the artificial viscosity, which reduces the artificial shearing 
stress. The explicit forms for the acceleration and the energy 
dissipation due to artificial viscosity for the i — eth SPH fluid 
element are: 

'^^ = -^m^nyV^Wij, (11) 
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where 11 is defined by (see, e.g., Lee and Ramirez-Ruiz 2002) 
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The function fi is defined by: 

f ^ |v-^l. 

|V--i;|,; + |Vxi;|, +77'c,//i,' 

and rj = 10^^. The sound speed and smoothing length of 
each element are denoted by Cj and hi respectively, and the 
factor r}' ~ 10~^ in the denominator prevents numerical di- 
vergences, afc = /3f, = 7/2 are constants of order unity and 7 
is the polytropic index from the equation of state. This form 
of the artificial viscosity supresses the shearing stresses when 
the compression in the fluid is low and the vorticity is high, 
|V-t>|^|Vxt;|, but remains in effect if the compression 
dominates in the flow |V-t;|^|Vxt;|. 

3.2. Initial Conditions and Applied Perturbations 

To construct a particular fluid torus, one needs only to spec- 
ify the equation of state, the distribution and absolute value 
of angular momentum, and the degree to which the potential 
well is filled (through the constant $0). We restrict ourselves 
here to systems in which t{r) = £0 = est. Thus the actual 
value of angular momentum fixes the position of the center of 
the torus, tq, as defined in 12.21 Numerically, this is done by 
means of a Monte Carlo technique, by distributing fluid ele- 
ments over the prescribed volume according to the analytical 
density profile, and subsequently relaxing them with artificial 
damping included in the equations of motion. The rotation 
law is strictly enforced during this procedure, and the inter- 
nal energy is fixed by the adiabatic relation so that one has 



complete control over the initial condition. We have verified 
that our initial conditions are indeed in equilibrium by evolv- 
ing them explicitly in time without any applied perturbations. 
No global evolution is observed during these computations. 

We have applied two fundamentally different types of 
perturbations to initial models: impulsive and periodic. In 
the first, an additional velocity field is imposed on the ini- 
tial condition at t = 0, and the system is evolved without 
additional perturbations for many (typically one hundred) dy- 
namical times. In the second, an additional, periodic term is 
added to the equations of motion to produce an acceleration. 
This can be applied only in the radial direction, or vertical, or 
both. 

Finally, one can hold £0 constant during a calculation, or 
vary it slowly (i.e., over many dynamical times). In the first 
case, the torus will remain essentially in the same radial re- 
gion, oscillating due to the applied perturbation. In the sec- 
ond, in addition to these oscillations it will either move in- 
ward or outward, depending on whether £0 decreases or in- 
creases. We have considered this possibility in view of the 
fact that gas within an accretion disk will, in reality, drift radi- 
ally as angular momentum is transported by viscous stresses. 

The temporal profile of the applied forcing can be varied. 
We have used single-frequency sinusoids as well as narrow 
pulse-like functions with repetition time Tg = l/i^s- This 
can be thought of as the rotation period of the central mass, 
which affects the accretion disk through some unspecified 
mechanism (e.g., the pulsar magnetic field, asymmetries in 
the quadrupole moment or other effects). 

4. Forced Oscillations 

We will hereafter restrict the presentation to the case of slen- 
der tori, where their radial and vertical extent, L, is small 
compared to the radial position of their center, i.e., L <C rg. 
The thin annulus can then be viewed as a small section of 
the entire accretion structure surrounding the compact ob- 
ject. The main difficulty with this picture, and in relating it 
to real systems, lies in the fact that the artifical zero-pressure 
boundaries obviously make wave propagation into the exte- 
rior of this torus impossible. Mode leaking is certainly an 
important issue (Fragile 2005), but for now we address only 
closed systems. The dynamical response of thick tori under 
global impulsive perturbations has been addresssed in a se- 
ries of papers by Rezzolla and collaborators (Zanotti et al. 
2003; Rezzolla et al. 2003a,b; Montero et al. 2004; Zanotti et 
al. 2005), while global and localized periodic perturbations 
have been considered by Rubio-Herrera & Lee (2005a,b). 

As a result of the induced forcing, the torus experiences 
small-amplitude oscillations, which can be measured in var- 
ious ways. One of these is to simply consider the position 
of its center, (ro,zo), defined, as above, as the location of 
maximum density, as a function of time. The corresponding 
time series are complicated, as they prominently display the 
perturbing frequency I's and the local epicyclic frequency for 
small radial and vertical oscillations, I'r and ly^ respectively. 
There are however, combination frequencies and cross-over 



4. 1 Radial and Vertical Mode Coupling 
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Fig. 3. Fourier transforms of the radial (solid) and vertical 
(dashed) oscillations of the center of a slender torus with 
constant specific angular momentum when perturbed peri- 
odically at frequency = 400 Hz in a pseudo-Newtonian 
potential. The local values of the radial and vertical epicyclic 
frequencies are j/^ = 500 Hz and I'z = 700 Hz, respectively. 
Even though the perturbation is purely radial, vertical oscil- 
lations are excited because of pressure coupling. The power 
was re-scaled along the vertical axis for illustrative purposes. 
In reality the power in vertical motions is much smaller than 
in the radial ones. 

phenomena, so that the whole behavior is best analyzed by 
performing Fourier transforms of ro{t) and zo{t). 

4.1. Radial and Vertical Mode Coupling 

For calculations where the angular momentum does not vary 
in time, since it is also constant in space, there are no mecha- 
nisms for its transport, and the perturbation does not produce 
any torques, the fluid must remain close to its equilibrium ra- 
dius tq. Figure |3] shows the Fourier transform of tq and zq 
for such a calculation, assuming the potential <E>kl, a cen- 
tral mass M = 1.4 Mq, tq = 6.1 ig and a purely radial, 
sinusoidal perturbation at frequency — 400 Hz. The cor- 
responding values of the local radial and vertical epicyclic 
frequencies are u,. = 500 Hz and I'z = 700 Hz. The power 
spectrum of radial motions clearly shows the perturbation 
frequency, i^g, and the radial epicyclic frequency, Vr- Like- 
wise, the power spectrum of vertical motions exhibits the 
vertical epicyclic frequency i/^. In this particular case, the 
value of To and is such that the difference between the two 
epicyclic frequencies is equal to half the spin frequency, i.e., 
~ Vr — Vs/2. There is thus clear evidence for mode cou- 
pling, through pressure, since the perturbation was initially 
only applied in the radial direction. Beats between the various 



frequencies are also present. In Figure|3] the power in vertical 
motions shows peaks at — Vg, v^, and Vg + Vj.. The very 
mode responsible for the coupling between the radial and ver- 
tical oscillations is weakly visible at — v,- = 200 Hz. 

A calculation with varying angular momentum (see ^ \3.2\ 
where the torus drifts inward is even more revealing. Start- 
ing at ro = 6.7rg, and terminating at tq — 5.35rg over two 
hundred initial orbital periods, the radial and vertical oscilla- 
tions of the center of the torus occur at a range of frequencies, 
covering values between the extremes in epiclycic frequen- 
cies Vr and Vz- This choice of parameters implies that, with 
Vg — 400 Hz, the calculation begins with ly^ — Vr < Vgjl and 
ends with ~ Vr > Vg/2. 

A power spectrum of the radial oscillations of tq dur- 
ing the a complete simulation using a pulse-like perturba- 
tion shows several features (see Figure]^. The perturbation 
at Vg is clearly visible, as are its three harmonics (it is not a 
pure sine wave, see Figure 3 in Lee, Abramowicz & Kluzniak 
2004). A broad, flat-topped peak is also apparent, covering 
the range 440Hz < Vr < 580Hz. These are simply the local 
values of Vr at the inner and outer radii, respectively. Har- 
monics of the radial epicyclic frequency, as well as a beat 
with Vg are also visible, at greatly reduced power. The corre- 
sponding power spectrum for the oscillations of zq shows a 
similar broad peak, for 590Hz < v^ < 880Hz, as expected 
(see Figure|3i. 

To explore the behavior of the oscillations in more de- 
tail, we have split the time series of t-q {t) and zq {t) into eight 
equal segments, each thus covering one eighth of the simula- 
tion. During this time, the torus is in a narrow radial window, 
and we may say that its position is nearly constant (except 
for the oscillations induced by the external forcing). We then 
extract the individual Fourier transforms for each segment. 
The results for the radial oscillations are shown in Figure |5j 
for the first and last time segments (they are all qualitatively 
similar). As expected, each one shows the perturbation and 
its harmonics, as well as a narrow peak centered at the value 
of Vr corresponding to the average position of the circle of 
maximum pressure during the time interval. The power spec- 
trum of the corresponding vertical oscillations is shown in 
Figurels] There is a strong peak at Vz, the amplitude of which 
varies greatly as the torus drifts radially. In Figure|6lwe show 
the power in this peak as a function of the average position of 
the center of the torus, rp. Two facts are immediately clear. 
First, the intensity of the coupling between radial and vertical 
modes is a function of the external perturbation. Second, this 
is a non-linear coupling, because the interaction is through 
a sub-harmonic of the pertrubation, Vg/2, and not with the 
perturbation itself or its higher harmonics. This fact points 
clearly to rich non-linear behavior in these systems, and, we 
believe, has been directly seen in accretion disks in LMXBs 
(Kluzniak et al. 2004; Lee, Abramowicz & Kluzniak 2004), 
in particular in SAXJ 1808.4-3658 (Wijnands et al. 2003). 

If the coupling between modes is due to some interac- 
tion with the external perturbation, and this excites a reso- 
nance, one would naively expect the consequences to display 
resonance-like behavior In particular, the resonant window 
should be narrow, and if the corresponding condition is not 
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met (e.g., above — Vr — Vs/2), then no coupling should 
occur. However, this is not what we observe in the calcula- 
tion we have just described. As the torus drifts inward, the 
power in vertical motions rises sharply as the resonance is 
approached, but remains high even after it has been passed. 
The corresponding range in frequencies for spans more 
than 100 Hz. This appears to indicate that once excited, the 
modes have locked (Kluzniak et al. 2004) and continue to 
be excited for some time. The curve in Figure |6l is reminis- 
cent in shape of the amplitudes of meridional oscillations for 
shghtly perturbed orbits in nearly geodesic motion, as com- 
puted by Abramowicz et al. (2003, see their Figure 2). In that 
case a generic coupling between radial and polar oscillators 
was assumed. In Neutron Star sources, the frequencies of the 
kHz QPOs drift over large ranges, sometimes hundreds of Hz. 
Mode locking following resonant excitation could in princi- 
ple be responsible for this behavior 

4.2, Additional Vertical Modes 

Inspection of Figures |3l and |5] reveals that the power spec- 
trum of vertical oscillations has a secondary peak at fre- 
quency Vu > Vz- The power in this oscillation is variable as 
the torus drifts radially, but it appears to be continously ex- 
cited. In some instances it seems to increase in power when 
Vz — Vr Vsl2, but its presence is not limited to the corre- 
sponding radial window. 

Thus it would seem again, in agreement with our previ- 
ous results, that a non-linear coupling between the radial and 
vertical modes transfers energy to the latter Further, a second 
oscillation at frequency ~ t^z + is excited within the 
same vertical motions, with power comparable to that at fre- 
quency Vz ■ Previously, when we examined peaks in the power 
spectrum at two different frequencies, they occurred in differ- 
ent variables (e.g., radial vs. vertical) and thus it was impos- 
sible to directly compare their amplitudes (as the strength of 
the real coupling is unknown, and the applied perturbation 
is arbitrary). Here we now show for the first time evidence 
for similar amplitudes in different modes of the same general 
motions. 

5. Conclusions and Directions for Future Work 

The oscillations of fluid confined by gravity and/or centrifu- 
gal forces can lead to rich behavior under a variety of as- 
trophysical scenarios (we have focused on QPOs in X-ray 
binaries, the observations of solar oscillations provide an- 
other, much closer example). Under quite general consid- 
erations we have shown that simple, acoustic modes which 
have their origin in free particle oscillations are modified by 
the hydrodynamics and can couple to one another in non- 
linear fashion. The locking of excited modes could explain 
the fact that the observed QPO frequencies drift over consid- 
erable ranges while still arising from resonant interactions. 
We believe that their signature is present in the observations 
of accretion flows in strong gravitational fields, and will al- 
low for the measurement of important parameters related to 
the central, compact object. 



Clearly, meaningful constraints will be derived only 
through much more detailed work. For example, it is not 
clear how such oscillations will actually modify the X-ray 
lightcurve (see, e.g., Schnittman 2005 for work along these 
lines). In addition, the strong gravitational field needs to be 
modeled in full General Relativity, with consideration of the 
self-gravity of the disk and detailed thermodynamics and ra- 
diation transport. 
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Fig. 4. Fourier transforms of the radial (top) and vertical (bot- 
tom) oscillations of the center of a radially drifting slender 
torus when perturbed periodically at frequency Us = 400 Hz 
in a pseudo-Newtonian potential. The outer (initial) and inner 
(final) central radii are ro = 6.7rg and ro = 5.35rg respec- 
tively. 
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Fig. 5. Fourier transforms of the radial (top) and vertical (bot- 
tom) oscillations of the center of a radially drifting slender 
torus when perturbed periodically at frequency Vs = 400 Hz 
in a pseudo-Newtonian potential. Only the first and last seg- 
ments of the calculation (1/8 of the total) are shown. 
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Fig. 6. Power in vertical motions at the local value of i^z as a 
function of the average central position of the torus, ro for the 
calculation with a radially drifting (inwards) perturbed torus. 
The calculation was divided into eight segments of equal du- 
ration. During each segment the torus drifted by 0. 16rg in the 
radial direction (shown by the horizontal bars). The strongest 
response occurs at ro ~ 6.15rg, where i/^ — = Note 
how the curve is not symmetrical with respect to the point of 
maximum power. See the text for details. 



